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ABSTRACT 



The magnetosphere of a pulsar is composed of relativistic plasmas stream- 
ing along the magnetic field lines and corotating with the pulsar. We study the 
intrinsic Faraday rotation in the pulsar magnetosphere by critically examining 
the wave modes and the variations of polarization properties for the circularly 
polarized natural modes under various assumptions about the magnetosphere 
plasma properties. Since it is difficult to describe analytically the Faraday ro- 
tation effect in such a plasma, we use numerical integrations to study the wave 
propagation effects in the corotating magnetosphere. Faraday rotation effect 
is identified among other propagation effects, such as wave mode coupling 
and the cyclotron absorption. In a highly symmetrical electron-positron pair 
plasma, the Faraday rotation effect is found to be negligible. Only for asym- 
metrical plasmas, such as the electron-ion streaming plasma, can the Faraday 
rotation effect become significant, and the Faraday rotation angle is found to 
be approximately proportional to A*^'^ instead of the usual A^-law. For such 
electrons-ion plasma of pulsar magnetosphere, the induced rotation measure 
becomes larger at higher frequencies, and should have opposite signs for the 
emissions from opposite magnetic poles. 
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1 INTRODUCTION 

Radio emissions from pulsars are generally highly polarized. When the radio waves propagate 
through the magneto-ionized interstellar medium (ISM), the linear polarization plane is 
rotated by an angle proportional to / Buedl ■ A^. This is the Faraday rotation (FR) effect. 
The rotation measure, RM = / Buedl, can be measured from the polarization angles at 
different frequencies (Hamilton & Lyne 1987; Han et al. 1999, 2006). 

Some questions need to be answered when pulsar rotation measures are used to inves- 
tigate the properties of the interstellar medium. Do the RMs of pulsars completely come 
from the ISM? Is there any intrinsic Faraday rotation from pulsar magnetosphere? If so, 
how large is it? Because the intrinsic polarization position angles of pulsar emission vary 
with the pulsar rotation phase, is the RM different for different phase bin? How to average 
the RMs by different weights of different phases? The RMs from ISM should not vary with 
rotation phase. However, recent observations by Ramachandran et al. (2004) and Noutsos et 
al. (2009) showed that for some pulsars the observed RMs do vary with the rotation phases, 
indicating that the RMs do not just come from the ISM. 

There are three possible reasons for the different phase-resolved RM: (1) the effect of 
incoherent superposition of quasi-orthogonal polarization modes of pulsar emission (Ra- 
machandran et al. 2004); (2) interstellar scattering (Noutsos et al. 2009, Karastergious 2009); 
(3) intrinsic Faraday rotation in pulsar magnetospheres. The first two possibilities have al- 
ready been discussed previously, while the third has not been addressed carefully. Although 
Kennett & Melrose (1998) discussed the generalized Faraday rotation effect in relativistic 
plasmas when the natural modes are linearly polarized, they did not consider the Faraday 
rotation for the circularly polarized natural modes in such a plasma. 

In this paper, we investigate the intrinsic Faraday rotation in pulsar magnetospheres 
when the natural wave modes are circularly polarized. A pulsar magnetosphere consists of 
highly magnetized relativistic streaming plasmas, different from the ISM. Thus, the simple 
relation for the Faraday rotation may not apply. In § 2, we give the theoretical description 
of the natural modes in the relativistic streaming plasma, and analyze the Faraday rotation 
effect. Because it is not easy to obtain analytic expression for the Faraday rotation effect for 
such a plasma, in § 3 we study this effect by numerically integration of the wave vectors in 
a self-consistent manner. Our conclusions are presented in § 4. 
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2 FARADAY ROTATION FOR CIRCULARLY POLARIZED NATURAL 
MODES IN PULSAR MAGNETOSPHERE 

In any magnetized medium, when the natural modes are circularly polarized and the mode 
evolution is adiabatic, Faraday rotation effect should rotate the linear polarization position 
angle (0pa) by 

/A 
—dr. (2.1) 

Here Ak = Anu/c is the wave number difference of the two circularly polarized natural 
modes, dr is the unit distance in photon ray ^. In our previous papers (Wang & Lai 2007; 
Wang, Lai & Han 2010, hereafter WLHIO), we have already derived the dielectric tensor and 
wave modes in pulsar magnetospheres composed of streaming, relativistic plasmas of various 
compositions and Lorentz factors. With the dielectric tensor, we can solve the Maxwell 
equation and obtain two eigenmodes, to be labeled as the plus "+" mode and minus "— " 
mode, and corresponding refractive indices n^-. We write the polarization vectors of the 
two modes as E± = E±t + E±^z in the xyz-iiamie, with z = B m. the xz-plane and 
k X B — — sin^Bif- The transverse part of the mode polarization vector is given by 

^±T = ^^^^|yT7^(^i^±,l), (2.2) 
where iK± — {Ex/Ey)±. From the Maxwell equation and dielectric tensor, we obtain 



X± = ^po, ± v//3p'oi + 1, (2.3) 
with the polarization parameter 

^ „i ~ fnSin'eBuj/il-u^r) ^2.4) 
2E^/i2,s(cos6'b - ^sSin6'B)' 

Here /i2,s, are factors in dielectric tensor [see eqs. (3.15) and (3.16) of WLHIO], 
the sum runs over each charged particle species "s" (electrons, positrons and/or ions). 
The dimensionless parameter Ur — ujc/[yuj (1 — ^cos^b)] stands for the ratio of cyclotron 
frequency cUc and the rest-frame frequency 70; (1 — P cos 9b), with cu the wave frequency, 7 

the Lorentz factor of streaming plasma, and 9b the k-B angle. Obviously, when |/3poi| ^ 1, 
the two eigenmodes are linearly polarized, while for |/3poi| ^ 1 the two modes are circularly 
polarized {K± — ±1). The refractive indices of the two modes can be written as 

n| = 1 + /ii ± K±J:Ji2,s{cos9b - C.sin^B). (2.5) 

^ In this paper, r stands for the distance from the neutron star center. Since the emission and propagation regions are far away 
from the NS surface, we do not distinguish the small difference of dr between the photon ray and radial direction. 
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Since n± is very close to the unity, the difference between the two refractive indices can be 
written as 

An ~ yr+^E,/i2,s(cos Ob - 6 sin Ob). (2.6) 

If the electrons and positrons in the plasma are purely symmetrical in density and veloc- 
ity, then the denominator of eq. (2.4) equals zero and ^poi = oo along the photon ray; this 
means that the two natural modes are always linearly polarized. However, in reality there 
always exists some asymmetry between electrons and positrons in the pulsar magnetosphere. 
For example, in the Goldreich- Julian magnetosphere model (Goldreich & Julian 1969), there 
exists a net charge density, called Goldreich- Julian density, A'^qj, which specifies the density 
difference between electrons and positrons. Ions (such as H+ or Fe ions) may also exist in 
the magnetosphere. Since the interaction between ions and photons is much weaker than 
that between electrons/positons and photons, the existence of ions enhances the asymmetry 
of the plasma. The Lorentz factors of electrons and positrons may also be different because 
of different electric field acceleration (Arons 1983; Kazbegi et al. 1991). If the asymmetry 
between electrons and positrons is sufficiently large, even though the natural modes in the 
inner magnetosphere of pulsars are linearly polarized, when the photons propagate to the 
outer magnetosphere, the natural modes could become elliptically or circularly polarized. 

In this paper, for simplicity we consider a cold pair plasma, in which electrons and 
positrons have the same Lorentz factors {'jp — — l) but different densities, AN — 
Np — Nf. ^ 0. The density difference could be Goldreich- Julian density, but we will not 
entirely restrict ourselves to this constraint. In this case all the equations are simplified. The 
polarization parameter /3poi becomes 

o ^ sin' ObT'' , . 

2(l-/3cos^B)(coseB-/3)AAr/Ar' 

with the plasma velocity ^ — l/Vl + 7^ and the total plasma density A^ = Np-\- N^. The 
difference of refractive indices can also be simplified to 

v^'^Ur AN cos6b — I3 



An- VI + /^poii + i27rad - ^^.^ 1-/3 costs' ^'"^^ 
Here v = Wp/cu' relates to the plasma frequency Wp, and the radiative damping parameter 

Trad = 7^ 1 2.9 

is important only near the cyclotron resonance and can be neglected away from the reso- 
nance. The cyclotron resonance radius rcyc can be written as [see eq. (4.49) of WLHIO] 



Tcyc 
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/R. = 1.8 X lQ^Bl{We-^'"r"'- (2.10) 

Here the magnetic field is described by B{r) ~ B^{R^/rY, is the radius of the neutron 
star (usually wc set 7?* = 10km), B^. = lO^^fi+12 G is the surface magnetic field, u = ug GHz 
is the wave frequency, and r is the distance from the neutron star center. Typically, 6b <^ 1, 
7^1 and 9bJ ^ 1, equations (2.7) and (2.8) can be simplified further: 

2Ur 



An.^l + /5^,-^^, (2.12) 

where we have neglected the radiative damping term. If » 1, the natural modes 

could become circularly polarized even before the cyclotron resonance (where u^- — 1), and 
Faraday rotation could affect the photon linear polarization angle. Before and after the 

cyclotron resonance (where Uj. = 1), An will change signs, while /3poi keeps the same sign. 
This means that the Faraday rotations to PA have opposite directions before and after 
cyclotron resonance. 



2.1 Faraday rotation effect for circularly polarized natural modes 

For a typical 9b (not so close to 0, 9b^ ^ 1), only when Ur < 9b^'^^N/N, the eigenmodes 
become circularly polarized [see eq. (2.11)]. We define the radius of circularization Tcir by 
\/3po\{r = rcn)\ = 1, so that the normal modes become circularly polarized when r > rdv 
According to eq. (2.11), the radius of circularization is 

n^/R. = 2.2 X 10'BlQP{;'iyg'^%'/'r\^N/Nr'/'. (2.13) 

For the typical parameters of a pulsar magnetosphere, B^ = lO^^G, P=ls, z/=lGHz, 9b ~ 
0.1, 7 = 100, we have Tcir/rLc ~ 0.1. If the mode evolution is adiabatic and refractive indices 
difference An is large enough, the linear polarization position angle will be affected by the 
Faraday rotation significantly. The position angle is rotated by 

where we have assumed f3poi ~ since the natural modes are circularly polarized. This 
integral cannot be evaluated analytically. The most important variable is B, since B oc r~^, 
and V (X B, Ur (X B. The other parameters in the equation do not change or change slowly 
with r. Thus, approximately we get 
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~ ^ ^TT^-y- / • 2.15 

Ur N 2C ' ^Wrcyc l-< V'^cyc/ 

Suppose x = r/rcyc, then ti^ — and the integration term of eq. (2.15) becomes 

' cir/ ' eye 

As discussed above, An changes signs when photon propagates across Tcyc (where x — 1). 
Then eq. (2.15) becomes 

A0PA ^ Q.l9riP{,'Bl%efr''\^N/N)F,v^''\ (2.17) 

where r] — N/Nqj is the plasma density parameter, Pig = P/ls. Equation (2.17) gives a 
simple approximate description to the final Faraday rotation angle. It is proportional to B^, 
77, AN/N, but inversely proportional to P, 7 and u. Notice that A(f)pfi^ oc Ug ^^"^ oc A^^^, quite 
different from the FR in ISM (for which A0pa oc A^). If the plasma is highly symmetrical 
{AN/N <^ 1), since 0pa oc AN/N, the Faraday rotation is negligible. However, if the plasma 
is highly asymmetrical, for example, AN/N ~ 1 (the case of pure electrons), rdi/rcyc could 
be less than 1, the Faraday rotation angle could be significant. For typical parameters of a 
pulsar magnetosphere, = 10^^ G, P = 1 s, 7 = 100, r] = 1000, AN/N ^1,63^ 0.1, the 
circularization radius can be rcir/Tcyc ~ 0.5, and the final rotation angle A(/)pA — 0.9ug 

Notice that in our derivation of equations (2.15) and (2.17), the k-B angle 63 is assumed 
to be a constant in the propagation. In reality, 9b also changes along the photon path. 
Additionally, given by eq. (2.16) is related to Tcir and Tcyc, both of which are determined 
by the plasma parameters [see eqs. (2.13) and (2.10)]. Thus, the Faraday rotation angle given 
in equations (2.15) and (2.17) is only a simple approximation. We will obtain the precise 
rotation angle by numerical integration of wave evolution equation in § 3. 



2.2 Evolution of mode amplitude and adiabatic condition 

In the xyz frame [defined with z — k, B — {— sin 9b, 0, cos 9b)], there are two wave modes: 
"+" mode and "— " mode. We introduce a mixing angle, 9m, via tan^^ = 1/{K+), so that 

tan2^^ = /3p-V (2.18) 

In the xyz frame, the transverse components of the mode eigenvectors are then 



E 



+ 



/ 








i cos 9„i 




' —i sin 9jn 1 




V 


sin 9m j 




^ cos 9m J 



(2.19) 



In the fixed observers' XYZ frame [defined hj Z = k, Q, in the XZ-plane and k x Cl ~ 
flsin^Y, the direction of B in this frame is {9b, <Pb)], they become 



E 



+ 



i cos 9m COS (pB — sin 9m sin b 

y i cos sin 05 + sin 9m cos 0^ 

The general wave amplitudes can be written as 
\ 

= A+E+ + A_E_. 



—i sin 9m cos 0b — cos sin 0b 
— i sin sin 0b + cos 9m cos 0b 



.(2.20) 



Ax 

Ay 



(2.21) 



Substitute this into the wave equation, we obtain the mode amplitude evolution equation: 



/ 



dr 



A+ 

A_ 



V 



(2.22) 



-Ak/2 + (f)'B sin 29m iO'm + 4>'b cos 29m 
-i9'^ + 0's cos 29m AA;/2 - 0'^ sin 2^„ 

where the superscript (') specifies d/dr, A/c = k+ — k_ = Anou/c, and a non-essential unity 
matrix has been subtracted. This equation generalizes the results in special cases (where 
only 9m or 0b varies) studied in Lai & Ho (2002, 2003) and van Adelsberg & Lai (2006), 
and therefore is useful for understanding the effect of mode coupling and Faraday rotation. 

When the natural modes are circularly polarized, which means 9m — 45°, eq. (2.22) can 
be simplified to 




-Ak/2 + cP'b i9'^ 



-i9L 



Ak/2-(f)'B 




(2.23) 



The adiabatic parameter is then defined as 



ad 



Ak/2 



91 



(2.24) 



Only when Fad ^ 1, the mode evolution is adiabatic, so that the two natural modes can 
propagate separately and the Faraday rotation can occur. For the non-adiabatic case of 
Tad "C 1, the two modes are coupled with each other, the electromagnetic waves do not 
interact with the medium and the wave polarization states keep unchange. 

For the wave propagation in the interstellar medium, the adiabatic condition is very 
easily satisfied. However, the pulsar magnetosphere is filled with relativistic streaming pair 
plasma, the condition is only satisfied in some regions. To evaluate the adiabatic condition 
Tad ^ 1 for r > Tcir (or (3po\ ~ 0), we use Ak = Ancu/c, with An given by eq. (2.12), and 

(I>'b = F^/rLC = 2.1 X 10-^°F^Pi-^ cm'^ (2.25) 



with 
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„ sin^ a cos C — sin a cos a sin C cos \I/ , , 

Fcp = ^ ^ 2- (2-26) 

1 — (cos a cos ( + sin ( sin a cos \E') 

Here a is the inclination angle between the magnetic axis and the rotation axis of a pulsar, 

C is the view angle from the rotation axis, and ^ = + is the pulsar rotation phase. 

Prom the definition of 9m, we have 

Here we define — l^-poi/ Ppo\ which is the scale length of ^poi variation; generally Lp r 
since ^poi oc r~^6'^. In the upper panels of Figure 1, we give some examples for the evolution 
of the quantities AA;/2, (p'g and O'^^^ along the photon ray for different plasma parameters. 
In the highly symmetrical plasma model (see Fig. la), the refractive indices difference An 
is small, and the evolution is non-adiabatic {Ak/2 ^ 0^, 9!^). Thus there is no Faraday 
rotation in this model. On the other hand, for the highly asymmetrical plasma model (see 
Fig. lb), the evolution remains adiabatic all the way along the photon path {Ak/2 » 0^, 
9'^). Therefore we expect that in this (asymmetric) plasma model, Faraday rotation will 
be generated. Only in some special cases, for example, if a photon propagates through the 
region where the magnetic fields change signs along the path very quickly, 0'^ could be larger 
than Ak/2, so that the non-adiabatic mode evolution could occur. This "quasi-transverse 
magnetic field" case was considered by Broderick & Blandford (2009) and Melrose (2010). 
However, the radio emission we considered in this paper is generated near the tangential 
direction of the magnetic field line, and magnetic fields in the outer regions do not change 
sign along the photon path, thus the "quasi-transverse" situation does not occur in pulsar 
magnetosphere. 



3 NUMERICAL RESULTS 

The analytic description above for the physical process of the Faraday rotation can only 
give qualitative results under simplified conditions. Here we use the numerical integrations 

to calculate how radio polarization evolves as the wave propagates through the pulsar mag- 
netosphere. This allows us to obtain the exact values for the rotated polarization angles 0pa 
and the rotation measure RM. 
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3.1 Single ray evolution 

It is generally accepted that pulsar radio emission originates from the open field line region 
at a few to tens of neutron star radii (e.g. Cordes 1978; Blaskiewicz et al. 1991; Kramer et 
al. 1997; Kijak & Gil 2003). We assume the emission height rem = 50i?^, and at the emission 
point, the photon is polarized in the k-B plane (i.e. the 0-mode in the case of curvature 
radiation) and propagates along the tangential direction of the local field hue. We neglect the 
emission cone of opening angle I/7. For a given emission height rem, pulsar initial rotation 
phase ^'i, the inclination angle a (i.e., the fx-fl angle), the view angle of line of sight ( (i.e., 
the k-Q, angle), the surface magnetic field 5*, and the plasma properties (plasma density 
parameter 77 = N/Nq;]. Lorentz factor of the streaming plasma 7), we can numerically 
calculate the dielectric tensor at each point along the photon ray [see equation (3.15) in 
WLHIO], and integrate the wave evolution equation [given by equation (2.10) in WLHIO] 
from the emission point to a very large radius, which we take the light cylinder radius ric, 
to determine the final polarization state of the photon. 

We first consider magnetospheres consisting of electrons and positrons, with the net 
charge density given by the Goldreich- Julian model, i.e., Np — N^. — Nqj. We assume that 
the total plasma density is N — Np + Ne — IOOOA^gj, which means the plasma is almost 
symmetric for electrons and positrons. Figure la shows an example of the photon polarization 
evolution along its trajectory in such a plasma. The wave natural modes are linearly polarized 
within the cyclotron absorption radius (~ 1650/2*). Before the cyclotron resonance, the 
dominated propagation effect is the wave mode coupling (around ~ 600R^), where circular 
polarization is generated. Although the natural modes become circularly polarized near 
r ~ 2600i?*, the refractive index difference is very small and the mode evolution is non- 
adiabatic too {Ak/2 < 0'^, 9!^); these make the Faraday rotation ineffective. Thus the PA 
remains unchanged outside the wave mode coupling area of ~ 600i?* (see the bottom panel 
of Fig. la). 

Figure lb shows the case where the magnetosphere plasma consists of electrons and ions, 
without positrons at all, N — — IOOOA^gj, — 0. The ions can be neglected for the 
emission and propagation effects since they are too heavy compared with the electrons. In 
this case, the adiabatic condition of Ak/2 > 0'^, O'^ is always satisfied (the upper panel of 
Fig. lb). The refractive index difference is also much larger than that for the symmetric case. 
In the region of r < 800i?*, the wave mode coupling dominates for the wave propagation. 
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Farther out, the natural wave modes becomes circularly polarized (rdr ~ 900-R*), Faraday 
rotation occurs and the direction of linear polarization varies (the bottom panel). Note that 
at the cyclotron resonance (rcyc ~ 1700i?*), Ak/2 changes sign, so the rotation measure 
changes sign too. The final Faraday rotation angle is about 20°. In Fig. Ic we consider the 
same case as in Fig. lb except for the emission coming from the other magnetic pole. 

3.2 Rotation measure in electron-ion plasma of pulsar magnetosphere 

As discussed above, the Faraday rotation effect is negligible in a highly symmetric electron- 
positron pair plasma of pulsar magnetosphere, while may be sigificant in an asymmetric 
plasma, e.g., electron-ion streams. In this subsection, we study the asymmetric case in more 
detail and present the rotation measure values for various plasma parameters. 

In Figure 2, we show some examples of the PA changes for different pulsar and magneto- 
sphere plasma properties, such as different electron density parameter 77 = N/Nqj (Fig. 2a), 
wave frequency u (Fig. 2b), surface magnetic field (Fig. 2c), pulsar period P (Fig. 2d), 
Lorentz factor of plasmas 7 (Fig. 2e) and the emission height rem (Fig- 2f). In general, the 
Faraday rotated PA increases with rj and S*, but decreases with u and P, and almost does 
not change with 7 and rem, roughly in agreement with eq. (2.17). 

Interpulses have been detected from a few tens of pulsars (e.g. Weltevrede & Johnston 
2008), which is 180" away from the main pulse and comes from the other magnetic pole 
of a pulsar. In this case, the magnetic fields along the photon trajectory have the opposite 
directions, and Faraday rotation should have the opposite sign. The PA evolution is also 
exactly opposite to that for the mainpulse case (see Fig. Ic). 

3.2.1 Faraday Rotation Measure 

For given pulsar and plasma properties, we can calculate the rotated PA for various fre- 
quencies. Note that some other propagation effects, such as wave mode coupling and quasi- 
tangential effect (see WLHIO), can also modify PA unless the modes evolution becomes 

non-adiabatic or the natural modes become circularly polarized. In general, the final PA 
value compared with the initial value at emission point is given by 

0PA = Mother + 0FR, (3.28) 

where 0fr is the Faraday rotation angle, and Mother is the rotation angle caused by other 
propagation effects in the magnetosphere. The Faraday rotation angles have the same mag- 
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nitude but different signs for the opposite magnetic poles. The other propagation effects, 
such as wave mode couphng, are the same for the opposite poles. Note that the cyclotron 
absorption effect will absorb opposite circular polarization but does not affect the linear po- 
larization angle. Therefore we have 0pA,opp — Mother — 0fr for the emission from the opposite 
magnetic pole. In general, we can compute 0pa and 0pA,opp using numerical integrations. To 
eliminate the PA rotation from other propagation effects, we simply use the 0pa value from 
the two poles to calculate Faraday rotation angle as 

0FR = (0PA - 0PA,opp)/2. (3.29) 

Figure 3 shows how the Faraday rotation angle varies with the wavelength. We find that 
is not proportional to in a pulsar magnetosphere (see the left panel of Fig. 3). For 
a set of given pulsar and plasma parameters, with surface magnetic field B^, = lO^^G, the 
best fit is (f)FR — 0.86 A°-^^ (see the right panel of Fig. 3). The index of 0.47 is different from 
the theoretical prediction of 1/3 in eq. (2.17), probably because of the simplifications of 
and 9b- In real observations, we always calculate the RM value by 

RM = (3.30) 

which varies with A in pulsar magnetosphere as shown in the middle plane of Figure 3. 
The higher the observed frequency is, the larger the RM is. Roughly, RM~ 0.2 A~^-^^ can 
describe RM dependence of A rather well. For the observation frequency of = 1.4GHz (or 
A ~ 21cm), the calculated RM^ 2 rad m~^. 

For a higher surface magnetic field, -B^, = 5 x lO^^G, the Faraday rotated angles are larger 
compared with those for — lO^^G (see Fig. 2c and Fig. 3). The relation between 0fr and 
A can be fitted to be 0fr = 2.28 A°-^^ The corresponding RM~ 0.64 A"^'^^. At i/ = 1.4GHz, 
RM~ 6 rad m^^ 

3.2.2 Phase resolved Faraday Rotation Measure 

Having studied the Faraday rotation angles and RMs for photons propagating through the 
magnetosphere with given initial pulsar rotation phases, we can now examine the phase 
resolved Faraday rotation effect. As shown in Fig. 2f, the emission height (much less than the 
light cylinder radius) does not significantly affect 0fr- Thus, for simplicity, we assume that 
all emissions of different rotation phases come from the same height, at rem = 50-R*. Figure 4 
shows examples of the Faraday rotation measure as a function of the pulsar rotation phase. 
To obtain RM, we integrate the photon polarization evolution across the magnetosphere. 
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subtract the PA variation caused by other propagation effects, and then get the Faraday 
rotated angle (pFR] we repeat this process with a shghtly different frequency, and finally 
obtain RM at each rotation phase. From Fig. 4, we see that the leading part of RMs are 
smaller than the trailing part. The maximum RM value appears at ~ 8°. The RM profiles 
change only shghtly for different impact angles. 

4 CONCLUSIONS 

Motivated by recent observations of phase-dependent Faraday rotation measures, we have 
studied the Faraday rotation effect of wave propagation in pulsar magnetospheres when the 
natural wave modes are circularly polarized. Pulsar magnetosphere is filled with relativistic 
streaming plasma, the Faraday rotation effect is different from that of the non-relativistic 
interstellar medium. We analyzed the wave modes and Faraday rotation effect considering 
the adiabatic evolution conditions. We used numerical integration of the polarization vector 
along the photon ray to incorporate all the propagation effects self-consistently within a 
single framework. We find that for highly symmetric pair plasma, with |A^p — A^e| ^ = 
Np + Ne, the magnetosphere Faraday rotation is negligible. For asymmetric plasmas (e.g., 
an electron-ion streams with ^ Nqj), the magnetosphere Faraday rotation effect may be 
significant. For such an electron-ion plasma of pulsar magnetosphere, the Faraday rotation 
angle 0fr is not proportional to in pulsar magnetospheres, but approximately to A°'^. 
The induced value of RM oc A~^-^, becomes larger for higher frequencies. The phase resolved 
RMs have smaller values for leading part and larger values for traihng part. There exists 
a maximum RM at the trailing half of the emission beam. The RM profiles do not change 
significantly for various impact angles. The induced rotation measures for the mainpulse and 
interpulse have the same magnitude but opposite sign. 

Our calculations in this paper have relied on several simplified assumptions. For example, 
we have assumed a cold streaming plasma with the Lorentz factor 7 and the density param- 
eter N/Nqj constant throughout the magnetosphere. Note that in reality, even when the 
Faraday rotation effect is very weak for a symmetric pair plasma, other propagation effects 
can modify the PA curve significantly (see Wang, Lai & Han 2010; Andrianov & Beskin 2010; 
Beskin & Philippov 2011). Overall, our result shows that an unambiguous identification of 
the magnetosphere Faraday rotation in the PA data may provide a valuable probe of the 
physical condition of pulsar magnetospheres. 
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Figure 1. Single photon evolutions across the pulsar magnetosphere with the plasma density parameter rj = N/Nqj = 1000 
and the Lorentz factor 7 = 100 for a pulsar with surface mag netic field B, = 10^^ q ^nd spin period P = 1 s. The magnetic 
inclination angle is a = 30° and the impact angle is x = C ~ = 5°. The calculations were made for the wave frequency 
u = 1 GHz, initial rotation phase VE"; = 0° , and emission height rem = 50-R* . Panel (a) is for the pair plasma with Np—Ng = Nqj . 



Panel (b) is the case for the magnetosphere plasma consisting of electrons and ions (without positrons), A'^ 
Panel (c) is almost the same as panel (b) but for the emission of the opposite magnetic pole. 
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Figure 2. The evolution of PA along the photon ray for various pulsar and magnetosphere plasma properties, including density 
Tj = N/Nqj (panel a), frequency (panel b), surface B field (panel c), rotation period (panel d), Lorentz factor (panel e) and 
initial emission height (panel f). In each panel, only one parameter is varied and labeled, while other parameters are the same 
as those used in Fig. lb. 
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Figure 3. The Faraday rotated angle produced in the magnetosphere as a function of the wavelength for = 10^^ G ("•") 
and B* = 5 X 10^^ G ("x"). The left panel shows (/>fr against \^ and the middle panel shows RM {=dif>p^/d\'^) against A. 
The right panel shows the possible relation between logio((/>FR) and logio(A), and the two solid lines are 0Fa = 0.86 A"'*'^ for 
S* = lO^^G and 0fr = 2.28X^'^^ for B* = 5 X lO^^G. The pulsar and plasma parameters are the same as used in Figure lb, 
except for the varying frequencies. 
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Figure 4. Phase resolved RM from pulsar magnetospheres for different impact angles of x = 5° (solid line), 2° (dashed line) 
and —5° (dotted line). The pulsar and plasma parameters for the calculation are the same as used in Figure lb, except for 
different initial rotation phase Stj and surface B field _B» = 5 X 10^^ G. 
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